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We study the effect of fluctuations on the ac conductivity of a layered superconductor both for 
c-axis and a&-plane electromagnetic wave polarizations. The fluctuation contributions of different 
physical nature and signs (paraconductivity, Maki-Thompson anomalous contribution, one-electron 
density of states renormalization) are found to be suppressed by the external field at different 
characterisitic frequencies (ojal ~ T — T c , u>mt ~ max{T — T c , T v x }, ojdos ~ min{T, r 1 } for the 
2D case). As a result the appearance of the nonmonotonic frequency dependence (pseudogap) in the 
infrared optical conductivity of HTS film is predicted. The effect has to be especially pronounced 
in the case of the electromagnetic field polarization along c-axis. 
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I. INTRODUCTION 
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Recently it was demonstratedi~EI that due to such specific properties of HTS as the extremely strong anisotropy and 
the small coherence length the fluctuation renoemalization of the one-electron density of states (or, in other words, the 
opening of the fluctuation pseudogap above T c )lteI can play an important role in the explanation of their normal phase 
properties near the superconducting transition. The depression of the one-electron density of states (DOS) on the 
Fermi level results in the decrease of related thermodynamic and transport properties which are proportional to that 
one. It is important that this effect has an opposite sign with respect to the traditionally discussed Aslamazov-Larkin 
QO \ (AL) and Maki-Thompson (MT) corrections though it remains less singular in T — T C E. 

The most explicit manifestation of the effect discussed can be observed in the c-axis transport properties of layered 
materials above T c , such as the fluctuation conductivity and the magnetoresistance, where paraconductivity andiJi 
Thompson contribution in the 2D regime turn out to be suppressed by the square of the interlayer trasmittanceErE 
Another interesting example of the evident appearance of the negative DOS contribution is the case of the NMR relax- 
ation rate for which the ALJike. contribution is rigorously forbidden for the singlet pairing and the MT contribution 
compets with the DOS oneli3il3. Strong pairbreaking weaks the former and the latter leads to the decrease of the 
$_i ' NMR relaxation rate at the edge of the transitiorila. 

In the present communication we attract the attention to another possibility to suppress the main AL and MT 
contributions by means of an high frequency electromagnetic field applied. Basing on the Drude origin of the normal 
part of the optical conductivity one can expect that noticeable changes in the electromagnetic wave reflectivity versus 
frequency take place for u> ~ r _1 . But it turns out that the frequency dependence of paraconductivity__jaiid Maki- 
Thompson anomalous contribution in the vicinity of T c appears long before, for u> ~ min {T - T^T-^Um. These 
positive contributions compet with the negative a (u>) which, as we will see below, weakly depends on frequency 
(only in the scale to ~ min{T, r -1 }). This competition results in the rapid decay of the dissipative processes already 
at frequencies of the order of u> ~ T — T c and in the appearance of the transparency window up to ui ~ min{T, t -1 }. 
The following high frequency behavior of Re [o"(w)] is mostly governed by cr n (w) which decreases, in agreement with 
Drude law, for u> > r _1 . 

Below we will study the ac fluctuation conductivity tensor for a layered superconductor taking into account all 
contributions and paying attention to the most interesting case of c-axis polarization of the field. Nevertheless, for 
completness, the afr-plane results of the old paper of Aslamazov and VarlamovEZI will be re-examined and discussed 
in application to the novel HTS layered systems. 



II. THE MODEL AND SOME DEFINITIONS 



We study the optical conductivity of layered superconductors which can be expressed by the retarded analitical 
continuation of the current-current correlator (electromagnetic response operator) (u>) : 



Re [(T a p{uj)] = -- 
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The problem of fluctuation conductivity of layered superconductors has been extensively studied and for the detailed 
description of the model and the diagrammatic presentation of the electromagnetic response operator tensor Q a p we 
refer tcO. So below we will list only the necessary definitions and results following the notations of that article. 

We assume the electron spectrum of the layered metal in the form 



f (p) = £ o(p) + J cos(p z s) - E F , 



(2) 



where £o(p) = P 2 /(2m), p = (p,p z ), P = (Px,Py) is a two-dimensional, intralayer wavevector, and J is an effective 
quasiparticlc nearest-neighbour interlayer hopping energy. The Fermi surface defined by the condition £(j>f) = is 
a corrugated cylinder, and Ep is the Fermi energy. This model spectrum is most appropriate for highly anisotropic 
layered materials, for which J/Ep <C 1. 

As it is well established now, the electron mean free path in HTS layered single crystals or epitaxial films turns out 
to be of the same order or several times larger than the coherence length £ a b(0), so the parameter T c r ~ 1 and the 
theory has to be constructed for an arbitrary impurity concentration. 

Intralayer quasiparticle scattering is included in the single quasiparticle normal state Green's function by means of 
the relaxation time r: 



G(p,u n ) = - 



1 



iu) n - £(p) ' 



(3) 



where u> n = w n [l +_L/(2|w„|t)]. 

The " Cooperon"Ej (triangle vertex accounting the interference in the impurity avaraging of the pairs of Green's £ue . c- . 
tions), which is necessary in the following calculations, may be expressed, for an arbitrary impurity concentration|13'tll 
as: 
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where 0(cc) is the Heaviside step function, and 
surface average with the spectrum ^ , we find: 

(K(p)-£(q-p)] 2 } = 



) denotes an average over the Fermi surface. Performing the Fermi 
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where vp = \pp\/m is the magnitude of the Fermi velocity parallel to the layers. In (Q), we have made the assumption 
tDQ 2 <C 1, which we suppose throughout this manuscript and its limits have been discussed ir£3. The fluctuation 
propagator L(q, oj^) in the vicinity of T c has the form : 




(6) 



where e = ln(T/T c ) » (T - T c )/T c for T - T c < T c , ip(x) is the digamma function, p = N(Q) = m/(2ns) is the 
single-spin quasiparticle normal density of states, and 
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is the positive constant which enters into the current expression in the phenomenological GL theory in two 
dimensionsEI 



III. PARACONDUCTIVITY 

Let us first examine the AL contribution (diagram 1 of Fig. |l|) to the ac fluctuation conductivity. The general 
expression for the appropriate contribution to the electromagnetic response operator as a function of the Matsubara 
frequencies of the external electromagnetic field uj v istJ: 

Q$(u>v) =2e 2 T]T / ^B a (q,uv,^)L(q,ug x 
J [ > 

(8) 



2 



The Green's functions block is defined as: 
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x G(p,w„ + „)G(p,w„)G(q- p,Wp_ n ). 
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In the vicinity of T c , for frequencies uj <C T, the leading singular contribution to the response operator Q a g 



arises from the fluctuation propagators in 
sufhccs to neglect its frequency dependence! 



rather than from the frequency dependences of the B a blocks, so it 
. This approximation leads to 
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Using these expressions in (||) followed by analitical continuation of the external Matsubara frequencies to the 
imaginary axis and integration over momenta, one can find the explicit expression for the imaginary part of the 
retarded electromagnetic response operatorB for real frequencies uj <C T: 
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where D 1 (z) = 2 In ^/i + y^ + r) , D 2 (z) = - y/z(z + r), ADi (z) = D x (z) - D x (e) , AL> 2 (z) = £> 2 (z) - £> 2 (e) and 

r = Arj J 2 /vp . The value r ~ £ 2 (0) / s 2 is the usual anisotropy parameter! characterizing the dimensional crossover from 
the 2D to the 3D regime in the thermodynamic fluctuation behavior at T c (except for a±, for which the crossover 
is from 0D to 3D at T c ). The expressions presented above solve the problem of the frequency dependence of the 
paraconductivity tensor in the general form for e <C 1 and uj ^ T for an arbitrary relation between e,r and u>, but 
they are too cumbersome. 

Let us concentrate on the most interesting case for the HTS analysis of 2D fluctuations where £, C (T) -€1 s (r <C e) 
and turns out to be suppressed by the necessity of the indo»G»dcnt tunneling of each electron participating in the 
fluctuation pairing from one CuOi layer to the neighbour oneOO. The approximation r <C e simplifies considerably 
the expressions ( |TT] ) and (|l^)u and it keeps, in the same time, their validity up to frequencies comparable to T c : 
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Let us mention two facts following from the expressions obtained. Firstly, the paraconductivity begins to decrease 
rapidly with the increase of frequency already for oj > T — T c (the critical exponents of this power decrease coincide 
with those ones in e-dependence of dc-conductivity tensor components: uu = 1 (2D fluctuations) and v±_ = 2 (OD 
fluctuations)). Secondly the assumption to neglect the w-dependence of the Green's functions blocks evidently breaks 
down at frequencies w > T and this dependence has the only effect to accelerate the decrease. 



IV. DENSITY OF STATES CONTRIBUTION 



As far as concern the DOS contribution to the electromagnetic response operator tensor the four main diagrams for 
it (other two of this kind are negligible in the case under consideration!^) are presented in Fig. [j] (2-5). The general 
expression for the DOS contribution to Q a p(u>) from diagram 2 is: 



(15) 
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)G 2 (p, w„)G(q - p, Lu^ n )G(p, 



The external frequency w„ enters the expression ( p"5| ) only by the Green's function G(p, u> n +u) and it is not involved in 
q integration. So, near T C1 even in the case of an arbitrary external frequency, we can choose the propagator frequency 
Wj, = 0. In the same way it can be treated the contribution from diagram 3 of Fig. |l|. 

The diagrams 2,4 and 3,5 of Fig. |l| are topologically equivalent and this fact would let think that they give equal 
contributions to <j(u>). Nevertheless the thoroughful analysis of the analitical continuation over the external frequency 
showsj that their contributions differ slightly and for the total DOS contribution to conductivity one can find: 
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Let us stress that, in contrast to ( pi] ) and (|l2|), this result has been found with the only assumption e <S 1, so it is 
valid for any frequency, any impurity concentration and any dimensionality of the fluctuation behavior. The function 
K [u>, T, t) can be easily used to fit experimental data. Nevertheless we present the asymptotics of the expression (|l7j ) 
for clean and dirty cases. In the dirty case 
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and in the clean case 
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for r" 1 < w < T 
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V. MAKI-THOMPSON CONTRIBUTION 

The total contribution of the MT-like diagrams to <j a p (oA-,has Jseen analysed i n0 in the case of zero frequency 
and the frequency dependence of a\\ (w) has been studied inE-3. Intll it was shown that actually the regular part of 
MT diagram can almost always be omitted. So we will not discuss it in this paper and we will concentrate on the 
anomalous MT contribution (Fig. ^ diagram 6). 

IntD it was demonstrated that in the case of quasi two-dimensional electron motion (^J) there is no formal necessity 
to introduce the pairbreaking time t v because the Maki-Thompson logarithmic divergency is automatically cut off 
due to the possible interlayer hopping. Nevertheless, all evidences show that the intrinsic pairbreaking in HTS 
is strong (at least one of its sources may be identified as thermal phonons) and an estimation of the appropriate 
t v ~ 2-=- 5 • 10~ 13 s is only several times larger than T c _1 . So we have to speak actually about the overdamped xegime 
for Maki-Thompson contribution and it doesn't manifest itself noticeably in the e dependence of conductivity^ (the 
major part of experimental results is explained in terms of AL or AL and DOS contributions). Anyway we are 
interested in MT contribution because of its frequency dependence which evidently determines another characteristic 
scale in addition to the previous three (T — T c , T, r" 1 ) we have introduced: wmt ~ tZ . 

We start, as usual, from the general expression of the anomalous MT contribution to the electromagnetic operator 
tensortll: 
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After the integration over p momentum and the summation over w„ in the range u) n G [— Wj,,0[ (anomalous part), 
one can find: 
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At this stage of calculations we artificially introduce the phase-breaking time in the "Cooperon" vertices. Carrying 
out the integration and separating the real and the imaginary parts we have: 
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where 7 = . In the case of two-dimensional overdamped regime (r <C e ^ 7) the expression ( p3| ) gives the 

following limits: 
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Let us remind that the expression (|23|) has been obtained without any limitation on frequency. Nevertheless we 
have done the assumption tDQ <C 1 through all the paper and it turns out that, while this condition doesn't 
restrict our results for the AL and the DOS contributions over the full range of frequency, temperature and impurity 
concentration, this is not so for the MT contribution. In fact, as it was shown inE£l, in the ultra-clean (or non-local) 
limit, when Tr > the assumption tDQ 2 <C 1 is violated for MT contribution and the results obtained are not 

valid there. Nevertheless one can see that this non-local situation can be realized in the clean case (Tr ^> 1) only 
and for temperatures in the range l/(Tr) 2 Ce<Cl. We suppose Tr ~ 1, as it is in the case of HTS, so we skip the 
discussion about the non-local limit. 



VI. DISCUSSION 



Let us start from the analysis of each fluctuation contribution separately and then we will discuss their interplay 
in Re [<7j_(u;)]. Because of the large number of parameters entering the expressions we restrict our consideration to 
the c-axis component of the conductivity tensor in the region of 2D fluctuations (above Lawrence-Doniach crossover 
temperature) . The in-plane component will be overviewed in the end of this section. 

The AL contribution describes the fluctuation condensate response to the electromagnetic field applied. The 
component of the current associated with it can be treated as the precursor phenomenon of the screening currents in 
the superconducting phase. Above T c the virtual Cooper pairs binding energy gives rise to a pseudogap of the order of 
T — T c , so it is not surprising that at higher frequencies the AL contribution decreases with the further increase of w. 
Actually w AL ~ T — T c is the only relevant scale for <r AL : its frequency dependence doesn't contain T, and r. The 
independence from the latter manifests mathematically the fact that elastic impurities do not represent obstacles for 
the motion of Cooper pairs. The interaction of the electromagnetic wave with the fluctuation Cooper pairs resembles, 
in some way, the anomalous skin-effect where its reflection is determined by the interaction with the free electron 
system. 

Another effect related with the formation of the fluctuation Cooper pairs, but on the self-intersecting trajectories 
(like the weak localization correction), is described by the MT anomalous contribution. Being the contribution related 
with the Cooper pairs electric charge transfer it doesn't depend on the elastic scattering time but it turns out to be 
extremely sensitive to the phase-breaking mechanisms. So two characteristic scales turn out to be relevant in the 
frequency dependence of that one: T - T c and r^ 1 . In the case of HTS, where 1 has to be estimated as at least 
0.1T C for temperatures up to 5 4- 10 K above T c , the MT contribution is overdamped, it is determined by the value of 
T0 and it almost does not depend on temperature. 

The density of states fluctuation renormalization gives quite different contribution to Re [ct(w)] with respect to 
those above. Physically it is related with the decrease of the one electron density due to the involvement of some 
number of electrons in the fluctuation Cooper pairing. At low frequencies (ui <C r _1 ) the lack of electron states on 
the Fermi level leads to an opposite effect in comparison with AL and MT contributions: Re [cr DOS (o;)] turns out 
to be negative and this means the increase of the surface impedance, or, in other words, the decrease of reflectance. 
Nevertheless, the electromagnetic field applied affects the electron distribution and at high frequencies u> ~ r _1 the 
DOS contribution changes its sign. It is interesting that DOS contribution, as one-electron effect, depends on the 
impurity scattering similiarly to the normal Drude conductivity. The decrease of Re [<7 DOS (k;)] starts already at 
frequencies u> ~ min{T, r -1 } which for HTS are much higher than T — T c and r^ 1 . 

The scenario of Refcr^ '] w-dependence with the most natural choice of parameters (r <C e < r^ 1 min{T, r -1 }) is 
presented in Fig. |2|. The positive AL and MT effects, in their w-dependence, are well pronounced at low frequencies 
on the background of the DOS contribution which remains in this region a negative constant. Then at u> ~ T — T c 



G 



the former decays and the Reoj_ remains negative up to us ~ min{T, r -1 }. The DOS correction changes its sign at 
u> ~ r _1 and then it rapidly decreases. The following high frequency behavior is governed by the Drude law. So 
one can see that the characteristic pseudogap-like behavior in the frequency dependence of the optical conductivity, 
mentioned in the title of this paper, takes place in the range u> G [T — T c , r -1 ]. The depth of the window increases 
logarithmically with e when T tends to T c , as shown in Fig. 0. In the case of a6-plane optical conductivity the two first 
positive contributions are not suppressed by the interlayer transmittance, they exceed considerably the negative DOS 
contribution in a wide range of frequencies and any pseudogap-like behavior is unlikely in a^(u>): the reflectivity has 
to be of the metal kind. The comparison between ( |l4| ) and ([l6]) shows that the compensation of the two contributions 
could only take place at loq ~ T / lne which is out of the range of validity of the AL contribution. nn 

Let us compare now the results of our calculations with the experiments available. The recent measurements£3>E2l 
of the c-axis reflectivity spectra, in the FIR region on Y BaiCu^O% single crystals, show the response of a poor metal 
with the additional contributions from IR active phonon modes (which we do not discuss here). With the decrease 
of temperature the c-axis optical conductivity decreases showing a transition from a Drude-like to a pseudogap-like 
behavior at uj ~ 180 cmT x . This gap grows deeper below 180 K without any abrupt change at the superconducting 
transition temperature T c = 80K . 

Such experimentally observed behavior of the optical conductivity is in qualitative agreement with our results. 
Really the suppression of the density of states due to the superconducting fluctuations in the vicinity of T c leads 
to the decrease of reflectivity in the range of frequencies up to uj ~ r _1 . The magnitude of this depression slowly 
(logarithmically) increases with the decrease of temperature but, evidently, at the edge of the transition it reaches some 
saturation because of the crossover to the 3D regime in fluctuations (where instead of ln(l/e) one has ln(l/r) — y/e, see 
(|l6|)). So no singularity is expected in the value of the minimum even in the first order of perturbation theory. Below 
T c the fluctuation behavior of (Ag) is mostly symmetrical to that one above T c (seeBa) and, with the further decrease 
of temperature, the fluctuation pseudogap minimum in the optical conductivity smoothly transforms itself into the 
real superconducting gap, which opens, in HTS, very sharply. Let us stress, that the independence of the pseudogap 
threshold from temperature appears naturally in our theory: it is determined by ujq ~ r _1 (see JlS|), (|l9) and Fig. 
|j). Comparing Fig. || e Fig. [|, one can easy see that the threshold doesn't move when the temperature changes 
but varies when the inverse of the scattering rate changes. As far as concern the numerical value of ujq, supposing 
T c t = 0.35 (which is the value for the scattering rate of the sample used in the experiment,undcr considerations, 
that is also in the experimental range of the inverse of the scattering rate TV 0.3 -j- 0.7c3'E3), one can see that the 
pseudogap threshold is of the order of 200 cm -1 , in agreement with the experimental data, even from a quantitative 
point of viewtBEa. 

Our theory is, strictly speaking, valid only in the vicinity of the critical temperature, where e« 1. Nevertheless 
the logarithmic dependence on e of the result obatained gives grounds to believe that qualititatively the theory can 
be valid up to e = ln(T/T c ) ~ 1, so for temperatures up to 200 K in the experiment discussed. So the theory is again 
in agreement with the experimental value of temperature 180 K until which the pseudogap is observable. 

In conclusion we have calculated the tensor of optical conductivity for layered superconductors. The pseudogap-like 
minimum of its c-axis component in a wide range of frequencies for temperatures in the vicinity of T c is found. Its 
origin is related with the fluctuation density of states renormalization which can be treated as the opening of the 
fluctuation pseudogap. These result are qualitatively, and in some aspects quantitatively, in agreement with the 
recent experiments on Y Ba2Cu^Og, samples. Further experiments, with more anisotropic samples like BSSCO single 
crystals, would be useful, because the effect should be more prononced. 
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FIG. 1. We report the leading Feynman diagrams contributing to the optical conductivity. The wavy lines indicate the 
fluctuation propagators; the thin solid lines with arrows stand for the impurity-averaged normal-state electron Green's functions; 
the shaded partial circles are the vertex corrections arising from impurities; the dashed curves with central crosses are additional 
impurity renormalizations. Diagram 1 is the Aslamazov-Larkin contribution, diagrams 2-5 are the corrections from the density 
of states renormalization and diagram 6 is the Maki- Thompson diagram. 
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FIG. 2. The plot shows the dependence of the real part of conductivity, normalized on the Drude normal conductivity, on 
ui/T, "R[a'(uj)] — Re [cr(w)] /cr n . The dashed line refers to the aft-plane component of the conductivity tensor whose Drude 
normal conductivity is = pe 2 TVp. The solid line refers to the c-axis component whose Drude normal conductivity is 
al = a\J 2 s 2 /vl. In this plot we have put Tr = 0.3, E F /T = 50, r = 0.01, e = 0.04, Tt v = 4. 
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FIG. 3. The behavior of the c-axis component of conductivity frequency dependence, for different values of temperature, is 
shown. The solid line refers to e = 0.04; the dashed line refers to e = 0.06; the dot-dashed line refers to e = 0.08. Tt — 0.2 for 
all the curves. The other parameters of this plot are the same used in Fig. 0. 
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FIG. 4. The plot shows the dependence of Re [a ±{uf) / a\] on lo/T for different values of Tr. The solid line refers to Tr — 
the dot-dashed line refers to Tr = 0.3; the dashed line refers to Tr = 0.2. The other parameters of this plot are the same 
in Fig. |. 
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